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We find two new classes of exact solutions to the Einstein-Maxwell system 
q-i of equations. The matter content satisfies a linear equation of state consistent 

$h | with quark matter; a particular form of one of the gravitational potentials is 

specified to generate solutions. The exact solutions can be written in terms of 
elementary functions, and these can be related to quark matter in the presence 
of an electromagnetic field. The first class of solutions generalises the Mak and 
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Harko model. The second class of solutions does not admit any singularities in 

<N ■ 

the matter and gravitational potentials at the centre. 
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1 Introduction 

The existence of quark stars in hydrostatic equilibrium was first suggested by Itoh 1 
in a seminal treatment. The study of strange stars consisting of quark matter has 
stimulated much interest in the last few decades since this could represent the most 
energetically favourable state of baryon matter. Matter consisting of u, d and s quarks 
may be the absolute ground state of matter at zero pressure and temperature as sug- 
gested by Bodmer-. It is expected that strange stars form during the collapse of the 
core of a massive star after a supernova explosion-. In regions of low temperatures and 
sufficiently high densities hadrons are crushed into quark matter with color supercon- 
ducting phases which occur in the dense cores of neutron stars as remarked by Alford 4 . 
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Consequently the core of a proto-neutron star or neutron star provides the appropriate 
environment for ordinary matter to convert to strange quark matter. Another possibil- 
ity is that a rapidly spinning dense star can accrete sufficient mass to undergo a phase 
transition to form a strange star. 

As the physics of ultrahigh densities for quark matter is not well understood, re- 
searchers restrict their attention to the phenomenological MIT bag model- - -. In the 
bag model, the strange matter equation of state has a simple linear form given by 

P = l -{p-AB) (1) 

where p is the energy density, p is the isotropic pressure and B is the bag constant. 
The quark confinement is determined by the vacuum pressure B in the bag model that 
equilibrates the pressure of quarks thereby stabilising the system. Studies of particular 
compact astronomical objects indicate that they could be strange stars with the quark 
matter- - —. A candidate for a strange star may have been observed using the deep 
Chandra LETG+HRC-S observations; Drake et al 1 - 7 suggested that the X-ray source 
RXJ1856. 5-3754 may be such an object. Sotani et al— have used observational data 
on gravitational waves to obtain the equation of state for quark matter. Harko and 
Cheng— considered collapsing strange matter in spherically symmetric fields. Yilmaz 
and Baysal— studied charged strange matter in rotating fields. The role of anisotropy, 
with the linear equation of state (pQ), was pursed by Mak and Harko 21 and Sharma and 
Maharaj— who demonstrated exact analytical solutions. 

In a recent treatment Mak and Harko 23 found a charged strange quark star under 
the assumption of spherical symmetry and the existence of a conformal Killing vector. 
In this paper we consider the Einstein-Maxwell system of equations with the linear 
equation of state (pQ) and apply them to strange stars. The existence of a conformal 
symmetry is not an assumption that we make. We demonstrate that exact analytical 
solutions to the field equations are possible that contain the Mak-Harko model—. In 
Section 2, we rewrite the Einstein-Maxwell field equations for the static spherically 
line element as an equivalent set of differential equations utilising a transformation due 
to Durgapal and Bannerji— . We then obtain a new set of differential equations with 
the assistance of the bag equation of state ([1]) for strange matter. On specifying an 
explicit form for one of the gravitational potentials, we obtain a first order differential 
equation in the remaining potential in Section 3. In Section 4 we find a new class of 
exact solutions to the Einstein-Maxwell system. The model of Mak and Harko^ 3 - is 
regained as a special case. In Section 5 we present a second class of exact solutions 
that satisfy the Einstein-Maxwell system. This category of solutions has the desirable 
feature of not admitting singularities at the centre. 
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2 Spherically symmetric spacetimes 



Since our intention is to study relativistic stellar objects it seems reasonable, on physical 
grounds, to assume that spacetime is static and spherically symmetric. This is clearly 
consistent with models utilised to study physical processes in compact objects. The 
metric of a static spherically symmetric spacetime in curvature coordinates reads 

ds * = _ e 2*(r) dt 2 + e 2X(r) dr 2 + f2( ^2 + ^2 (2) 

where v{r) and A(r) are the two arbitrary functions. For charged perfect fluids the 
Einstein-Maxwell system of field equations are obtained as 

(3a) 
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for the line element (T5]). We are utilising units in which the coupling constant 8ttG/c 4 = 
1 and the speed of light c = 1. The energy density p and the pressure p are measured 
relative to the comoving fluid 4-velocity u a = e~ u 5^ and primes means a derivative 
with respect to the radial coordinate r. In the system ([3]) the quantities E and a are 
the electric field intensity and the proper charge density respectively. 
It is convenient at this point to introduce the transformation 

A 2 y\x) = e 2v{r \ Z(x) = e" 2A(r) , x = Cr 2 (4) 

where A and C are arbitrary constants. With this transformation, the system ([3]) has 
the equivalent form 

2 •• , o rr 2 • ^ — — - -E 2 ^ 



4Zx^y + 2Zx J y + ^a; - Z + 1 - —J y = (5c) 

^ = — (xE + E) 2 (5d) 

where dots denote differentiation with respect to the variable x. Note that equation 
fl5c|) is the condition of pressure isotropy. We can replace the system of field equations 
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dSJ), including the bag equation of state ([!]), by the system 
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The system of equations ([6]) governs the gravitational behaviour of a charged quark 
star. 



3 Integration procedure 

We describe one possible integration procedure that leads to an exact solution of the 
Einstein-Maxwell system ([6]). Note that other procedures are possible; our approach 
has the advantage of producing a first order equation that has solution in terms of 
elementary functions. We observe from (I6al) that p and p are related. Therefore in 
the system ([6]) there are five independent variables (Z, y, p or p, E, a) and only 
four independent equations. We have freedom to choose only one of the quantities 
involved. In our approach we specify y(x) on physical grounds. A number of choices 
for the gravitational potential y(x) are possible; clearly we should choose a form that 
is likely to lead to a physically reasonable solution. To make the above set of equations 
tractable, we choose the metric function in the particular form 

y(x) = (a + x m ) n (7) 

where a, m and n are constants. The form chosen ensures that the metric function y 
is continuous and well-behaved in the interior of the star for the wide range of values 
of parameters m and n. The function y yields a finite value at the centre of the star. 
This is a very desirable feature for the model on physical grounds. It is interesting to 
observe that many of the solutions found previously do not share this feature. 
Substitution of ([7]) into fl6d|) leads to the first order equation 

^ a 2 + 2a[l + mn(2m + l)]x m + [2mn(2mn + 1) + l]x 2m 
+ 2x(a + x m ) [a + (1 + mn)x m ] 

(1-f x){a + x m ) 
2x[a + (1 + mn)x m ] 
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This first order equation is linear so that it can be integrated in principle. The compli- 
cated rational coefficient of Z can be simplified using partial fractions and we obtain 



Z + 



1 2m(n - l)^" 1 
2x (a + x m ) 



m[4(l + mn) — 3n]x 
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Note that we have essentially reduced the solution of the field equations (JBj) to inte- 
grating (|S]). Once the potential Z in ([8]) is found the remaining relevant quantities p, p 
and E then follow from f l6al) . ( I6bl) and f l6cl) respectively. It is possible to find exact 
solutions to the Einstein-Maxwell field equations with the linear equation of state for 
different values of m and n in . We illustrate this with two simple examples in terms 
of elementary functions. Other exact solutions are possible but the form of the solution 
becomes more complicated and could involve special functions. 



4 Generalised Mak-Harko model 



An exact solution of ([6]) can be found with m = 1/2 and n = 1. In this case ((Tj) gives 
the first metric function 

y(x) = (a + s/x) 

Equation (jSJ) becomes 
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which can be integrated to give the second metric function 
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Hence we can generate the exact analytical model 
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which is a solution of ([6]). For simplicity we have set 

3C(6a 2 + 10a^+3x) 



2^(a + ^x){2a + 3y/x) 2 



The exact model ( fTUl) satisfies the Einstein-Maxwell system ([6]). Note that if we set 
a = 0, then the system (flOl) becomes 



e^^V, e 2 ^-^, P=-k+B,P=-k~B, E 2 = ± (li; 

1 — ±)r z 2r z or z 3r z 



where we have set D 2 = A 2 C . The particular solution (11 ip was found by Mak and 
Harko 23 under the assumption of spherical symmetry and the existence of a one- 
parameter group of conformal motions. It is interesting to observe that on substituting 
the values a = and B = for the constants, the solution fflUj) becomes identical to 
that obtained by Misner and Zapolsky 25 . The physical features of the solutions fill I) 
were studied by Mak and Harko 23 and shown to be consistent with the interior of a 
quark star with charged material. This corresponds to a single stable quark configura- 
tion with radius R = 9.46 Km and mass M = 2.86M ; these figures are consistent with 
values obtained using numerical methods by other researchers^ - —. Consequently our 
more general class of solutions is likely to produce charged quark models consistent 
with stellar evolution and observational data. We comment that our new class of so- 
lutions fflOl) has a singularity in the charge density and mass density at the centre; the 
Mak and Harkc-2 3 - model also shares this feature. The singularity in the charge density 
and mass density is physically acceptable since the total charge and mass remain fi- 
nite. However our gravitational potentials e 2v and e 2A remain finite at the centre which 
contrasts with the singularities in the metric functions of Mak and Harko when x = 0. 

5 Nonsingular quark model 

Another exact solution of (J6]) can be found with m = 1 and n = 2. For this case ((Tj) 
gives the first metric function 

y{x) = (a + x) 2 . 

Equation (jSJ) becomes 



Z + 



1 2 3 

7T- + + 7T- 
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- ^§x)(a + x) 
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which can be integrated to give the second metric function 

9(35a 3 + 35a 2 x + 21a:r 2 + 5x 3 ) - ^a;(105a 3 + 189a 2 a; + 135ax 2 + 35a; 3 ) 
~ 315(a + x) 2 (a + 3a;) ^ ' 
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Therefore we can find the exact analytical model 
e 2v = A 2 {a + xf (13a) 

315(a + x) 2 (a + 3x) 



e 2A 



9(35a 3 + 35a 2 x + 21ax 2 + 5x 3 ) - ^x(105a 3 + 189a 2 x + 135ax 2 + 35x 3 ) 

(13b) 
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P = ^ ) + 105(a + x) 3 (a + 3x) 2 (13c) 
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Again for simplicity we have set 

6C(70a 4 + 217a 3 a; + 159a 2 x 2 + 75ax 3 + 15x 4 ) 



9{x) 



35(a + x) 3 (a + 3x) 2 

The exact model ([TBI satisfies the Einstein-Maxwell system (jSj) and constitutes 
a new family of analytical solutions for a quark star with charged material. The 
gravitational potentials e 2u and e 2A in ([TBI) have the advantage of having a simple 
analytic form, and they are written in terms of polynomials and rational functions. 
Consequently the matter variables and the electric field intensity have a simple analytic 
representation. The function e 2u is continuous and well behaved in the interior and 
finite at the centre x = 0. The function e 2A is well behaved and has a constant value at 
the centre x = 0. The energy density p is positive throughout the interior, regular at 
the centre with value po — 2(— + B). The pressure p is regular at the centre with value 
Po = 2( 2 ^- — j) = |(p — AB). The electric field intensity E is continuous in the interior 
and vanishes at the centre. Hence the matter variables and gravitational potentials 
comply with usual conditions for a stellar source. The finiteness of e 2u , e 2X , p, p and 
E at the origin x = is a very welcome feature which is absent in the previous class of 
solutions. Consequently the exact solutions ([TBI) are likely to produce charged quark 
stars with physically acceptable interiors. A recent attempt in this direction is the 
strange star model of Jotania and Tikekar 29 admitting compact configurations with 
mass to size ratio consistent with strange matter. 
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6 Conclusion 



We have generated a new category of exact solutions to the Einstein-Maxwell system 
of equations. The linear equation of state (pQ) was imposed which is relevant in the 
description of quark stars. Two classes of exact solutions were identified by specifying 
the form of one of the gravitational potentials. The first class comprises the generalised 
Mak-Harko model as the earlier solution of Mak and Harko 23 is regained as a special 
case. The second class comprises the nonsingular quark model which has finite values 
for both the matter and metric variables at the stellar centre. The method in this 
paper depends crucially on the choice ([7]) for the metric potential Z which leads to 
solutions of the condition of pressure isotropy. In future work it would be desirable to 
seek physically reasonable solutions, with new forms of Z, which are consistent with the 
relationship ([1]). We also intend to study more closely the physical features of quark 
stars, from the exact solutions generated, and relate these to specific astronomical 
objects following the treatment of Jotania and Tikekar— . 

In conclusion, we make two further points that are of significance in the modelling 
process. Firstly, in the general solutions (TTUT) and (|T3|) . when studying models of 
charged spheres, we should consider only those values of parameters for which the 
energy density p, the pressure p and the electric field intensity E 2 are positive. Clearly 
a wide range of charged spheres, with nonsingular potentials and matter variables, are 
possible for relevant choices of a, B and C. Secondly, the interior metric ([2]) must 
match to the Reissner-Nordstrom exterior spacetime 



, 2M Q 2 \ , , / 2M Q 2 \ -1 , ., ., 
ds 2 = - 1 + -^)dt 2 + (1 + -— r- ) dr 2 + r 2 (dd 2 + sin 



2 , Jl( Ja2 i „;„2 

2 I '" 

across the boundary of the star. This generate relationships between the mass M, the 
charge Q, and the constants a, B and C. Consequently the continuity of the metric 
coefficients across the boundary of the star is easily achieved. There are sufficient free 
parameters available to satisfy the necessary conditions that may arise from a partic- 
ular physical model under consideration. 
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